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THEORY OF WINGS IN H0N3TATI0HARY FLOW* 
By A. I. HekraBov 


As long as the velocities of airplanes were relatively small 
aerodynamics was a science of steady motion In an Incompressible 
fluid. It was found that neglecting the compressibility of the air 
an! considering only steady motion led to results which, for small 
velocities of the airplanes, oould be applied In practice. The ' 
considerable Increase, particularly in reoent times, of the speeds 
of airplanes has forced attention on unsteady motion In Incompress- 
ible and compressible fluids. 

The most Important effeot on the development of nonstationary 
aerodynamics was that exerted by the phenomenon of flutter which 
appears in certain oases at large flight speeds. Failure of the 
airplane structure due to flutter have occurred In all countries. 
Survivors of flutter oatastrophles and casual observers have 
reported that on increasing the speed of fie airplajie a vibration ■ 
of the wing or tall surfaces began after which the wing or tall 
broke suddenly, the speed and force of the rupture being similar to 
those of an explosion. The phenomenon of flutter may be readily 
observed also In a wind tunnel In which is placed a thin vooden 
flat strip. If through a window In the wall of the tunnel this 
strip Is moved crosswise to the tunnel on Increasing the velocity 
of the flow, the strip at first begins to flutter and then at a 
oertain flow velocity its amplitude of vibration Increases so 
rapidly that the inexperienced or nonalert experimenter rarely 
succeeds In saving the Btrip from failure by removing it from the 
flow In time. The phenomenon of flutter Is more energetic than that 
of resonance. To obtain flutter there Is not required any outside 
vibratory motions for the wing or tall surfaces. The cause of the 
flutter of the wing or tail surfaces is evidently a purely aero- 
dynamic one but to understand this phenomenon qualitatively may be 
possible even from a consideration of the motion along a straight 
line of a material point. Let us assume that the material point 
with mass m moves along the axis .0 X under the action of farces 
the projection X of which on the 0* axis Is equal to 

X » - mk 2 ! - m2u 
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Then the equation of notion of the material point vlll he 

The Integral of this equation. If p Is small. Is of the form 
i «= ae"^ sin ( ^k 2 - p 2 t + c) 

where a and £ are arbitrary constants. We see that If the 
quantity p Is positive the material point will he In harmonlo 
vibratory notion with damping. We shall now assume that 11 Itself 
depends on a certain parameter w so that p ■ p(w), where for 
w » Wcr p = 0 and far w>w Qr p Is negative, that Is, equal 
to - X (where X > 0) . Then far w a w cr the above equation assumes 
the farm 

x » a Bln (kt + £ ) 

that Is, we shall have the usual harmonlo vibration of a point, but 
for v>w cr we obtain 

x = ae** sin ( •/hfi - X 2 t + e) 

that is, the amplitude of the vibration begins to Increase exponen- 
tially, Thus, whereas In the case of resonance the amplitude of 
vibration has, as Is known, the form at, in the case of flutter 
it has the form ae . Table I shows the increase of the amplitude 
with time In resonance and in flutter for different values of the 
quantity X . 
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20,086a 
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4a 

1,492a 

7,389a 

54,598a 

5 

5a 

1,649a 

12,183a 

148,413a 


We note that whatever the value of the positive number X the 
ratio e^/t approaches Infinity with increasing t. For the air- 
plane the parameter w is the velooity of flight where the velocity 
w a Wqp, the value at which the coefficient p becomes zero, Is 
oalled the critical flutter velocity. 
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It Is not difficult to understand the physical cause of the 
phenomenon of flutter. Let us assume the wing of the airplane fcr 
any reason whether, for example, a gust of air, acquires a very 
small vibration. The elastic forces evidently are always damping 
forces which tend to extinqulsh the vibrations of the wing or tall. 
The aerodynamic forces however can be both damping forces for certain 
types of vibrations and for cob which assist the vibrations. As long 
as the velocity of the airplane is small the elastic damping foroes 
are predominant over the aerodynamic forces. With inoreaee in the 
flight speed the elastic foroes renaln unchanged In magnitude but 
the aerodynamic forces Increase, as is known proportionally to the 
square of the velocity of flight. There must therefore be an 
instant at which the aerodynamic foroes, which assist the vibrations, 
become predominant over the elastlo foroes and this Is the Instant 
at which flutter begins. The problem of the designer is to construct 
the wing such that the critical velooity of the flutter is higher 
than the velocities at which the airplane will fly and whloh is 
sufficiently strong to withstand the vibratory action of the aero- 
dynamic forces. From this it is clear that even before constructing 
the airplane it ia necessary to be able to determine theoretically 
the critical velocity of fluttor of the a* - *-1006 or^y from its 
design and by changing the dc.» .gn to rc-.l bo tho critical flutter 
velooity. 

It is not difficult with the aid of simple examples to show 
that the aerodynamic forces due to the deformations of the wing or 
tail surfaces depend on tho magnitude of these deformations and tho 
rate of change of these deformations and this explains the necessity 
for the presence of terms in x and dx/dt in the differential 
equation of the motion of the material point, which was taken above 
as an example to illustrate the phenomenon of flutter. Let us assume 
first that the wing of the airplane received only a twist deforma- 
tion so that the angle of attack of the wing increased somewhat. 

Since the aerodynamic force supporting the wing Increases with 
Increase (up to a certain limit) of the angle of attack, as also the 
torsional aerodynamic moment, we conclude that the additional aero- 
dynamic force obtained on twisting the wing depends on the magnitude 
of this wing deformation and increases together with this magnitude. 
We assume next that the wing of the airplane reoolved only a bending 
deformation occurring at velooity v (fig. l) . As long as there wbb 
no bending deformation the angle of attack of the wing was a = ZMAB 
where w *= AB is the flight speed of the airplane. With the flight 
velooity w = AB and the velooity of bending deformation v => AC 
we shall have a total velooity equal to AD and the angle of attack 
becomes equal to ZMAD = a - Aa where tan Aa ■=■ v/w or on account 
of the smallness of Aa approximately Aa = v/w. To the deorease 
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of the angle of attach a by the amount Aa & v/v oorreeponds a 
decrease In the aerodynamic forces depending as, we have seen, on 
the velocity v of the deformation and not on the magnitude of 
this deformation. For downward bending of the wing with velocity v 
the angle of attach a changes Into the angle a + v/w as a result 
of vhioh the aerodynamic forces Increase. From what has been said 
It follows that the aerodynamic forces obtained In the deformation 
of the wing depend on the velocity v of this deformation and are 
damping forces so that they oppose the development of the deforma- 
tion. 


In computing the critical flutter velocity it is neoessary to 
deal both with the elastic properties of the wing or tail surfaoes 
and with the magnitudes tending to produce flutter. As regards the 
elastic forces, experience shows that up to 60 percent of the 
breaking loads the deformation of the wing Is proportional to the 
force, that Is, the wing behaves like a beam. For this reason the 
matter Is simple as regards the elastic forces of the wing and the 
results of the theory of elastic beams can be applied. 

As regards the aerodynamics of the fluttering wing the matter 
Is not however so simple. There Is as yet no definite clearness 
In this field of Investigation. It Is evident that the aerodynamics 
of nonsteady flow by the essential nature of the phenomena must In 
addition to flutter havo application In the theory of flight acro- 
batics, in studying the effects of gusts of air. In the theory of 
flapping wings and In the flight of birds, in which Leonardo da Vinci 
had already been Interested, etc. As a result of this It can he 
understood why the great Interest which- at the present time appears 
throughout the world in the study of the aerodynamics of nonsteady 
flow. In Germany, for example, about 10 years ago a special division 
of the Prandtl Institute was set apart to devote itself entirely to 
nonstationary aerodynamics, the Institute being headed by Kussner. 
This institute has already Issued a series of Important Investigation 
results. The American Glen-Martin company applies nonstet ionary 
aerodynamics to compute the airfoil flutter. In the dosing remarks 
of one of its publications the company points out how Important it 
is at the present time to direct all efforts toward the development 
of the aerodynamics of nonsteady flow. 

The aerodynamics of steady flow was first applied to the non- 
steady motions as for example, to wing flutter, that is, it was 
assumed that at eaoh instant at every position of the wing the flow 
about the latter was such as though the wing were In this position 
an infinitely long time. As it turned out, the application of the 
aerodynamics of steady flow to nonsteady conditions of motion gives 
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In a number of oases results that are useful for application, and 
this has permitted the use -of the aerodynamics of steady flow up to 
the present time. With further Investigations, however, It became 
clear that In a number of other oases the application of the aero- 
dynamics of steady flow to nonsteady motions gives wrong results. 

The air forces acting on various types of wing profiles under 
steady flow conditions wore studied by Joukovaky, Kutta, grand tl, 
Chaplygin and others In Investigations many of which even now are of 
fundamental Importance. The study of the forces seeing on wings 
under nonsteady flow condition: presents far the general case of the 
nonsteady motion as oxprosscd by Proud tl, "a problem of transcendental 
difficulty". But it is precisely these particular cases of nonsteady 
flow which present the greatest Interest for aeronautics that have 
all been found amenable to investigation. 

ThuB was explained that the motion of an airplane in a circle 
with constant angular velocity, as approximately occurs In acrobatic 
flights on entering a dive, proceeds with constant circulation as 
a result of which there nay not be a vortex vn\e behlnu the vlng. 

The first one to study this type of motion was S. A. Chaplygin in 
1326 and after him Independently Glauert In 1929, completing his 
Investigations with the addition of a vortex wake behind the wing. 

On the other hand in the case of the phenomenon of flutter of the 
wings and tail surfaces of an airplane it was found sufficient to 
investigate the phenonemon only for small vibrations and this enables 
the investigation cf problems with variable circulation. 

The first investigator who laid the foundation for the present 
day study of the aerodynamic forces acting on a vibrating wing was 
Bimbaum (1924), a student of Prandtl. He introduced t^e important 
conoepts of free and bound vortices and vortex wake, Kussner developed 
the initial ideas of Bimbaum and thus advanced the problem of the 
airfoil in nonsteady motion. Beginning in 1935 them appeared a 
number of Hussion papers the authors of which were mainly Keldish, 
Sedov, and .Lavrentyev who indicated ways in which we can apply the 
theory of the functions of a oomplex variable to the problems of 
the nonstationary dynamics of an airfoil. In 1941 N. S. Kochin gave 
a strict solution of the problem of the v J o-rations of a wing of 
ciroular plan form, making only the esev.i option of linearity of the 
equations. The character of the solution in oloBed form as given 
by Kochin appears the only one at present for a wing of finite span. 

In 1936 Karman and Sears in the United States concerned themselves 
with problems of nonsteady flow and gave a new method of approach to 
the problem, thereby founding the American school of investigators. 

In following this school we study as it were the instantaneous photo- 
graph of the aerodynamic phenomena as they ooour. The methods of the 
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German school on the contrary deal as It were vlth the kimsmato- 
graphlo picture of the phenomena since they deal vlth the development 
of the phenomena In time. In 1938 Cicala In Italy made the first 
attempt to extend the theory of the lifting vortex of Prandtl -jjo the 
case of a wing of finite span In a nonsteady flow. Recently Kussner, 
starting from the method of acceleration potential introduced by 
Prandtl returned to the problem of Cicala. In the same year there 
appeared In Italy the work of Posslo who attempted to solve the 
problem of the nonsteady aerodynamics of a profile for the case of 
a compressible fluid. In his investigations Posslo makes use of the 
acceleration potential of Prandtl. Thus at the present time there 
are three independent schools occupying themselves with the problems 
of the aerodynamics of unsteady flow: the American, Russian, and 

Gorman. 


As in the usual aerodynamics of steady motion the problem of 
the nonsteady motion of the Infinite span airfoil is most advanced 
while there hardly exists any ocmplete theory of the unsteady motion 
of a wing of finite span of arbitrary plan form even for the 
Incompressible fluid although for a wing of circular plan form there 
exists the completed solution of Kochin. 

By the nonsteady motion of a wing there is understood In all 
these theories the following. It is assumed that the wing has a 
certain finite velocity w constant in magnitude and direction; in 
addition to the velocity w the wing possesses, as an absolutely 
rigid body, arbitrary infinitely small deviations from this velocity. 
Finally, the wing need not be absolutely rigid but may have infinitely 
small deformations varying with time. 

In order to render the problems under consideration linear it 
is assumed that the wings are infinitely thin, the profile in the 
absence of deformations coinciding with a linear segment. The 
angles of attack of the wing are assumed infinitesimal, in the 
fundamental equations or In the boundary conditions there being kept 
only the terms of first older smallness . Only with these restric- 
tions imposed does the problem of the unsteady motion of a wing 
become solvable at the present time in general form. What Is 
required to find in the solution of these problems is the magnitude 
of the foroes and moments which act on the wing in its unsteady 
motion. 

It was found by Joukovsky that the existence of a lift force 
on a profile is explained by the presence around the profile of 
. circulation on the magnitude of which the lifting force depends. If 
this profile moves with constant velocity the circulation is constant 
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but as a result of the nonuniformity of the motion of the profile 
the circulation around the profile must be variable since the lift ' 
force for nonuniform' motion of the profile varies. - Since the total 
circulation by the lavs of hydrodynamics cannot change there must 
appear behind the profile In nonuniform motion a vortex wake 
consisting of the continuously shed vortices. As has been said, 
Bimbaum was the first to turn ■ attention to the existence of the 
vortex wake. The vortices forming the vake are not displaced behind 
the moving profile but remain motionless In the fluid at rest. If we 
neglect their small displacement at right angles to the direction of 
the velocity of the airfoil. 

In all these theories there are considered principally periodic 
motions since It can be shown that an aperlodlo phenomenon may be 
represented as the sum of an Infinite number of Infinitely small 
periodic phenomena with different frequencies. Let (i) be the 
angular frequency of vibration of the airfoil, c half the profile 
chord and w the forward velocity of motion of the airfoil. The 
nondlmenslonal number k = (Oc/w Is called the Strouhal number which 
Is constantly applied In these theories. 

It Is first of all necessary to solve the problem to what 
extent the reotlllnear vortex wake can be considered stable since 
all theories assume Its rectlllnearlty. The Investigation of this 
problem Is rather a delicate one but In any case It may be carried 
out for a single particular case where the profile Is not absolutely 
rigid. In this case it is founi that If the vertical velocity on 
the profile is small the amplitudes of the vibrations of the vortex 
wake about its rectilinear form will likewise be small, the value of 
the velocity of the wake perpendicular to Its length being propagated 
as a wave along the wake with the velocity w = coc/k and on the 
vortex wake there will be waves with these velocities and 
frequency 0 ). 

According to the views of the American school there exists at 
each Instant on the profile a quasi-statlonary vortex intensity 7 q 
which is the vortex Intensity on the profile if at a given instant 
all the vortex wake were rolled into a point at infinity. Since for 
each Instant of time there will be a value of this Intensity although 
In all cases we Imagine the vortex wake rolled up Into a point at 
Infinity the Intensity Is called a quasi -stationary vortex Intensity 
In contrast to the stationary vortex Intensity which oocurs in the 
absence of a vortex wake and whose magnitude Is constant In time. 
Since the vortex wake at all times exlstB directly behind the 
trailing edge of the moving airfoil up to infinity, the total vortex 
intensity on the airfoil is equal to 70 + 71 wbere 71 is due to 
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tbs vortex wako. Barman shoved that knowing the funotions jtq and 
7 ]_ there may he- found the vortex intensity in the wake . From this 
it is possible to find the momentum and angular momentum of a fluid 
due to the entire system of vortices on the profile and in the 
vortex wake. Knowing the momentum and angular momentum of the fluid 
surrounding the profile it is possible to find by the theory of 
Euler the pressure force L and the moment M of the pressure 
forces on the airfoil. Karman ■ shoved that the force L and the 
moment M . may be represented as the sum of three components : 

L * Lq + Li + Ig, M * Mq + Mi + Mg 


The force Lq is the quasi -stationary force of Joukovsky to 
which corresponds the moment Mq. The magnitudes 1 q and Mo nay 
exist alone only in the limiting case when the airfoil moves with- 
out acceleration and there is no vortex wake behind the airfoil. 

The force and moment M]_ are due to the associated mass of the 
profile. To obtain them it is necessary to assume that the profile 
moves nonuniformly but without tho presence of circulation about it 
so that there 1 b no vortex wake behind the profile. Finally, the 
force Lg and the moment Mg are due to the vortex wake, where 
knowing the expressions for these it is not difficult to prove the 
following theorem due to Glauert : The resultant of the pressure 
foroes on an infinitely thin straight profile that is due to the 
vortex wake is applied at the forward foous of the profile, that is, 
at one quarter chord from the leading edge. 


In the theory of Barman and in the other theories there is used 
the function C(k) of the Strouhal number k introduced by the 
American investigator Thoodorsen. The function C(k) of The odor sen 
is of the form 


C(k) - F(k) + iG(k) 


Hi ( 2) (k) 

Hj/ 2 )(k) + iHo^ 2 ) 00 


where the functions E*/ 2 ) (k) and Hg^) 00 aro the Hankel func- 
tions satisfying the equation of Bessel and ha- lng the form 




00 
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A tat 1 b of! funotions of Theodoreen far the real values of the 
Strouhal number k has been set up. whore C(O) a 1. Although there 
have- been objections against, dividing, the £arae....L. and moment M 
In which It was pointed out that this splitting up was artificial 
these objections do not appear to us as well founded slnoo such a 
splitting up appears also possible in the -work of the German and 
Russian school and has moreover an entirely concrete physical basis. 

Proceeding to the consideration of the Russian school on the 
theory of the airfoil- In a nonsteady flow It Is first of all necessary 
to dwell on the introduction of funotions of a oompjLex variable into . 
this theory. As is known, the description of Bteady motion of an 
Incompressible fluid by functions -of a complex variable is entirely 
natural and has boen productive of important results. Joukovsky and 
Chaplygin conducted their vork on airfoil theory in a steady flow 
using the theory of functions of a complex variable. Thus the work 
of the Russian school on airfoil theory in nonsteady flow may be 
considered as tho natural continuation and extension of the work of 
Joukovsky and Chaplygin. In this field however wo have considerable 
difficulties since the functions describing the unsteady phenomena 
depend not only on the complex variable z = x + iy but also on a 
roal variable, namely, tho time t. It was possible to overcame 
these rather Important difficulties and obtain an extension of the 
theorem of Joukovsky and Chaplygin on the lift farce and the moment 
of the pressure forcos for an arbitrarily deformed airfoil in a non- 
lncompresslblo fluid for tho qases of both fixed and moving 
coordinate axes. .Expanding the derivative the complex potential in 
powers of the complex variable z in the most general form, assuming 
that the coefficients of this expansion depend on the time t and 
introducing a certain auxiliary function g(z) it was possible 
with the aid of .the extended theorems of Joukovsky and Chaplygin to 
obtain expressions for the forces and moments acting on the airfoil 
in unsteady motion. The obtained expressions for the force L and 
moment M as in tho theory of tho American school were obtained by 
following the somo principle, that is, - splitting up into three 
components : 


L = Jjq + L]_ + M — Mq + M]_ + Mg 

The expressions far the quantities Lg and Mg obtained by 
the Russian school wore found to be identical with the expressions 
obtained for Idle analogous quantities by the Amorlcan school. As 
regards the quantities Lq, and Mq, obta±ned_ by the 

Russian and Amorlcan schools there is however an essontlal difference. 
In tho formulas of the American school those quantities are expressod 
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In terms of the quasi -stationary Intensity 70 whereas In the 
formulas of the Russian sohool these components are expressed In 
terms of the component v y (x, t) of the Telocity of the fluid near 
the airfoil perpondicular to the forward veloolty w of the air- 
foil. The transformation of the formulas of one school Into those 
of the other has as yet not been effeotod. 

The work of N. E. Kochin on the aerodynamics of a vibrating 
wing of circular plan form could not of course be carried out In 
terms of functions of a complex variable since It is not a two- 
dimensional but a three-dimensional problem. This work appeared as 
a result of an earlier work by Kochin on the steady motion of such 
a wing. It Is to be noted In these works by Kochin among others 
that Kochin was able to construct a function possessing definite 
properties and singularities by which It was possible to express 
the solution in closed form without making use of series. The 
form of the circle was taken far the reason that only for It, has 
It as yet been possible to find the above mentioned special func- 
tion. All restrictions and assumptions which Kochin Introduced In 
his Investigations reduoe to those by which the problems are made 
linear. Tor this reason the solutions of Kochin both for the steady 
and nonsteady flow are the most accurate of the existing solutions. 

A comparison of the results obtained for concrete problems on steady 
flow by the theory of Kochin and according to general approximate 
theories shows that the difference may amount to several tens per- 
cent of the magnitude to be determined. Tram this it Is clear how 
Important It Is to have exact solutions of the type of the solution 
of Kochin. It Is desirable therefore that by the same method It be 
possible to solve the problem of the elliptic wing which resembles 
mare nearly an actual -wing than does the circular wing but at the 
same time It Is Impossible to close our eyes to the fact that such 
a solution presents extremely great difficulties and requires great 
mathematical aptitude. 

Although behind the profile In unsteady motion there should 
as a general rule exist a vortex wake It may be proven that even 
for periodic motion of a profile there may be particular oases of 
Its motion where there will exist no vortex wake behind the profile. 
This will be the case for example when the projection of the velocity 
at the profile at Its rear focus, that Is, at one quarter chord from 
its trailing edge in a direction at right angles to forward veloolty w 
of the profile Is equal to zero, and also in the case where the 
absolutely rigid airfoil moves over a cirole with constant angular 
velocity. In the first case the circulation about the profile must 
be equal to zero and In the second case It must In general be constant. 
In 1926 Chaplygin considered a number of profiles with oonstant 
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circulation about them In circular motion as absolutely rigid bodies. 
In 1935 the work of Cbaplygln was continued by Sedov by the extension 
^ oy, rangg of, profiles and bya .detailed- analysis -of -the analytic 
* structure- of the expressions for the associated mass and associated 
moment of inertia. As has already .been said above, in 1929 Glauert 
arrived at the same results Independently of the work of Chaplygin 
but then Glauert supplemented his work by considering elliptic 
profiles with vortex wakes behind them making uses of the results 
already available at that time to take aooount of the effect of the 
vortex wake . 

For studying the oases of motion of an airfoil with constant 
circulation about it both Chaplygin and Glauert make use of a moving 
system of coordinate axes rigidly fixed to the airfoil and obtain 
the equations of hydrodynamics with the Integral of Lagrange 
determining the pressure of the fluid at any point. From thlB, after 
rather complicated computations, it is found possible to find the 
expressions for the projections on the axes of coordinates of the 
resultant pressure and the resultant moment of the pressure foroes 
exerted by the fluid on the airfoil. In particular, Chaplygin 
applied the results obtained by him to the elliptic airfoil and to 
the rectangular airfoil moving over a circle with constant angular 
velooity. This problem is an example of unsteady motion differing 
essentially from rectilinear and not differing only by infinitely 
small deviations from the latter. Thus, assuming the circulation 
oonstant it is found possible to reject the condition of infinitely 
small disturbances from a uniform rectilinear motion. 

We now proceed with the presentation of the results of the 
German sohool of investigators. If, around a sufficiently thin 
rectangular airfoil ( -c, +c) there exists a constant circula- 
tion, the vortex intensity T) on the airfoil is oonstant, that 
is. Independent of time, and by the usual formulas of Helmholtz it 
is possible to compute the velooity at any point of the fluid due 
to the vortex intensity distribution on the airfoil. The case will 
be otherwise however if this intensity T| depends on the time t 
since in the first plaoe there will exist behind the airfoil a vortex 
wake which will have an effect on the value of the required, velooity 
and in the second place the free vortices forming the vortex wake 
behind the airfoil will exist on the airfoil itself so that to the 
intensity T) of the vortices on the airfoil must be. added the 
intensity e of the free vortices at the profile. The formation 
of the free vortices at the airfoil is explained by the fact that 
the. increase At| during each infinitesimal o'lement of time' of the 
intensity tj on the airfoil must be accompanied by the formation 
during the same time Interval of free vortices with the inten- 
sity - Atj. The free vorticeB forming about eaoh point of the 
airfoil, which remain fixed in space at the point at which they 
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are formed , will continuously maintain their intensity until, thei 
airfoil passes by- them. When the airfoil In Its ' displacement will 
coincide with them the latter, remaining stationarily, will oontinue 
the vortex wake behind the airfoil. Thus also for unsteady motion 
the projections of the velocity at any point of the fluid can be 
computed by the formulas of Helmholtz but their expressions will 
consist of three components depending respectively on the vortex 
intensity in the wake, • the vortex intensity e of the free 
vortices at the airfoil and the vortex intensity c' of the free 
vortices in the vortex wake. It Is not difficult to find the 
analytical expressions that connect the functions e and £' with 
the funotlon tj. Hence the projection v y of the velocity of a 
fluid on a direction perpendicular to the vector w can be expressed 
In terms of the function i), the expression consisting of three 
components containing the funotlon T] under the integral sign. Such 
expression for the projection Vy . of the velocity of the fluid was 
first obtained by Birnbaum. Assuming that the motion of the airfoil 
and of the fluid are periodic with frequency 0), that Is, Introducing 
for the functions depending on the time the exponential factor e* ) ' b 
we can divide all the obtained expressions for the projection Vy by 
the quantity e 1 * 0 *. It 1 b evident that this expression for the 
projection v y can be applied to the fluid particles In the imnedlate 
neighborhood of the airfoil. Sinoe the projection Vy for the 
particles at the airfoil can also be determined from the klnematlo 
considerations while the vortex Intensity T) an the airfoil Is an 
unknown function, the obtained equation may be looked upon as an 
integral equation for the function T). The latter is called the 
Integral equation of Birnbaum. The solution of this equation Is 
required since the lifting force, that is, the resultant force L of 
the pressures on the airfoil and the resultant moment H of the 
pressure forces exerted by the fluid on the airfoil are expressed in 
terms of the function tj, The solution of the Integral equation of 
Birnbaum is obtained In two forms. In the first place there are 
found the relations by which the coefficients of the series for the 
required function T) are expressed in terms of the coefficients of 
the series of the given funotlon Vy. In the second place there Is 
obtained a solution of the equation of Birnbaum In closed form, that 
Is, a formula Is found expressing the funotlon T) In terms of the 
funotlon Vy, 

In order to give an Idea of the analytical complexities which 
are here encountered we present both the equation of Birnbaum and its 
solution In dosed form. We assume that the airfoil (-o, +c) Is 
placed along the Ox axis with Its center at the point 0. We set 

x o c cos 0 or x = o cos p; tj = f (0) e 1 ^, v y = w(0) e 1 ^ 
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The Integral equation of Blrnbaum Is of the form 


&.(,.) .[* - yj [* 

Jo COS (1 - COB 0 iJo 


Q-lk QOS 9 s j n Pdg 

cos \i - cos 9 


F* f (<p) e 1 ^ 008 *P sin cpdcp + f 

Jo Jo 


“ e-Ik oh'A ah Ad* 
cob (i - oh X 


J ' 11 f(<p) e 1 ^ 000 ^ sin qpdcp^- 

whore k Is the Strouhal number, k = toc/w. The solution of this 
equation In closed form Is 

t(e ) = | -J"* ^[- COB £ + CM (1 + COS £ )J 1 9 + 

+ f .in e lg [ y. ^ ■.fl±|i]\v(£) d£ + | 

2 L l - cos ( e - e)J J * 


r* [V(B) - *(£)] -syeu- 

j 0 u coo 6 - cos e 

vhere C(k) Is the function of Theodorsen. 

Having the expressions for the force L and the moment M In 
terms of the function t] and making use of the solution of the 
equation of Blrnbaum It was found possible to obtain the expressions 
for the force L and the moment M In terms of the coefficients 
of the series of the function v y . The magnitudes L and M can 
be decomposed into three components following the ideas of the 
American school. It may be noted that the lift force L = Lq + Lq + 1«2 
depends only on the three first coefficients and the moment 
M = Mq + Mi + M2 depends on the coefficients A q, Aj., Ag, A3, A4, 

Ag . , . of the expansion of the function v y In a trigonometric 
series. In particular, for an absolutely rigid rectangular airfoil 
in an unsteady motion it was found that the lift force of the fluid 
is reduced to three definite forces applied at the forward focal 
point, at the center and at the rear focal point of the airfoil and 
to a pair of forces with definite moment . 
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If the airfoil 1b a rectangular absolutely rigid segment the 
mathematical analysis In determining the resistance oan be brought 
to completion. It Is found that the resistance force Is due In this 
case not only to the formation of the Barman vortex street behind 
the airfoil but to the formation behind the airfoil of a vortex 
vake. Whereas the vortex street of Barman always leads to a 
resistance, the force due to the formation of a vortex wake may be 
positive as well as negative, that Is, a thrust or a drag. In the 
case of a thrust It Is important of course to know the efficiency. 

The efficiency for the translational motion alone and the 

efficiency rj r for the rotational motion alone can in any case be 
determined. In these computations It is necessary to take into 
account the suotlon force obtained as a result of the pressure of 
great reduction In pressure at the leading edge of the rectangular 
wing. In figure 2 are shown the values of these coefficients ritr 
(the top curve) and ri r (the lower curve) as functions of the 
Strouhal number k. 

We see that the coefficient utr which Is important for small 
values of the Strouhal number k drops to half Its value on 
increasing the number k; the coefficient is always positive 

for all valueB of k, that is, for only translational fluctuations 
the airfoil alwayB possesses a thrust force. The coefficient rj r 
however Increases with increase in the Strouhal number. At k = 0.954 
approximately the coefficient n r becomes zero. Thus, for k < 0.954 
the rotating ving has only a resistance force but for k > 0.954 It 
has a thrust force, that ie, the latter is obtained only for rapid 
rotational motion of the airfoil. Although these conclusions were 
obtained for an airfoil of infinite span it is useful to compare 
them with the flight data of birds whose wings perform predominantly 
translational fluctuations in flight (table 2). 



a 

H 

M 

too" 1 

m/cm 

k 

Stork 

0,500 

2,03 

0,123 

12,6 

n 

0,134 

Pigeon 

0,075 

0,67 

0,045 

31,4 

25 

0,056 

Book 

0,130 

0,95 

0,068 

19,0 

12 

0,108 

Seagull 

0,126 

0,96 

0,066 

13,9 

12 

0,076 


In the above table S is the area of the wings, L their 
span, c the mean half chord equal to S/2L, to the angular 
frequency, v the flight speed, k the Strouhal number. We see 
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that "birds fly at small Strouhal numbers for which their motion 
has a considerable efficiency. Recently V. Golubev has worked out 
a method of competing the thrust force of an. airfoil in translational 
oscillatory motion on the basis of the existence of two v o rtex 
streets whose origin is at the extreme positions of the translational 
fluctuations of the airfoil and extending behind the airfoil. The 
direction of rotation of the vortices in these vortex streets is 
opposite to the rotation of the Karman vortices. Thus the vortices 
of Golubev push the fluid behind the airfoil and thereby move the 
airfoil forward. In this case infinitely small oscillations of the 
airfoil are not required. 

The methods developed by the German sOhool may also be applied 
to a oampoBite airfoil consisting of an absolutely rigid airfoil, 
aileron and trimmer. It is then necessary to compute a large 
number of magnitudes. These computations vere carried out twice by 
two different’ methods by Kussner in Germany and Theordorsen in 
America. From the agreement between the obtained results it follows 
that these results can be reliably used. 

The results obtained for periodic phenomena may, on applying 
the Fourier integral, be extended to phenomena which are not 
periodic in time. The mathematical analysis is here rather complicated, 
a fundamental part being played by the function k-^a) given by 



where the path of integration consists of the entire real axis with 
the origin bypassed below over an infinitely small ssmicirole. In 
order to give an example of how rapidly these aperiodic phenomena 
are developed we assume that the velocity v y at the airfoil at 
the Instant t = 0 has received a certain Instantaneous increase. 
Then the force L and the moment M likewise change their values 
disoontinnously hut only by half their final value. If the chard 
of the airfoil is 4 motors and the veloolty w = 100 meters per second 
then at the end of 2 seconds the force L is 99 percent of its final 
increment. The total Increase of the force L and the moment M 
will occur after an infinitely long time, that is, to attain the 
remaining 1 peroent an infinite time is required. 

In the same as the phenomena which are aperiodic with respect 
to time we can consider the phenomena which are aperiodic with 
respect to space. It is' possible far example to' consider" the case 
where the airfoil during its motion enters a rising current of air. 

We can obtain the solution by considering the system of infinitely 


■16 


HACA TM Ho. 1154 


many solutions with different Strou&al numbers. In tbls theory the 
function 

* 3 (o) - [*,(«) - IJltO] 0(.) 4. 

Jv-* 

■where Jq(z) and J^(z) are Bessel functions and C(z) the 
Theodorsen function, plays an Important part. In order to obtain 
an Idea as to how rapidly the airfoil reacts to the effect of the 
rising air current we may point out that after the airfoil has 
moved a distance equal to fifty times Its chord In the rising 
ourrent the lnorease In the lift force attains 99 percent of Its 
total value . 

It has been pointed out above that the extension of the problem 
of the unsteady motion of an airfoil to the case of an Incompressible 
fluid was made by the Italian Fosslo. In addition to general 
restrictions by which the problems are rendered linear Fosslo 
introduced Into the problem a further restriction namely that the 
velocity of sound a In the gas Is always a constant quantity whereas 
for adiabatic processes the velocity of sound depends on the square of 
the velocity of the gas and a constant value for it is only approximate. 
To solve the problem Fosslo makes use of the potential acceleration 
of Prandtl, the partial derivatives with respect to the coordinates of 
which are equal to the corresponding projections of the accelerations. 

It may be shown that for the assumed restrictions the acceleration 
potential satisfies the same seoond order partial differential 
equation as the velocity potential. At Infinity the acceleration 
potential must be equal to zero and on the airfoil have a certain 
discontinuity. Erom this It follows that the acceleration potential 
may be represented as a potential of dipoles distributed along the 
airfoil with axes perpendicular to the latter. We may observe that 
this dipole does not In Its dimensions correspond to an actual dipole 
since the word Is here used to denote a particular solution of the 
equation for the acceleration potential. The analytical expression 
for the latter will contain under the Integral sign the intensity 
of these dipoles. From this analytical expression It Is possible to 
find the accelerations by simple differentiation and then the 
velocities of the compressible fluid, the transition from accelera- 
tions to velocities, under the assumptions and restrictions made, 
being effected through the Integration of the linear equations. We 
can then find the projection of the velooity of the gas on the normal 
to the vector w In the Immediate neighborhood of the airfoil and 
the expression for this projection will contain under the Integral- 
sign the unknown intensity of the dipoles. On the other hand this 
projection may be determined directly on the basis of kinematic 
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considerations . Equating the two expressions we obtain an Integral 
equation for determining the unknown dipole intensity. Since the 
lifting force . L and total moment M are expressed in terms of 
this intensity the solution of this Integral equation furnishes the 
solution of the problem. Fosslo succeeded in obtaining this 
Integral equation but neither he nor anyone else has succeeded as 
yet in solving this equation. The equation has the following form 



where p, is the Mach number equal to w/a, p ■ to/w, the 
function is a Henkel function and the symbol |ti( represents 

the absolute value of tj. From this equation it is required to 
determine the intensity q(x') of the dipoles on the airfoil for 
given value of v y (x) of the projections of the velocity on the 
airfoil. From the mere appearance of the equation it is clear how 
difficult this problem must be. 

The theory of the wing of finite span in an unsteady flow was 
treated by two different methods, the method of Kussner, using the 
acceleration potential and the method of Cicala introducing a speolal 
system of vortices consisting of one infinite rectilinear lifting 
vortex and three systems of infinitely many rectilinear and 
reotangular vortices. After complicated explanations. Cicala, for 
determining the circulation r(y^) at any perpendicular seotion yi 
of the wing, arrives at the following integro-diff erential equation 


▼z(n) 


where 


n(fc) r(yi) 

2no 


■ j * ,) 




yi - y 


2s being the span of the wing, 2c its chord, k * coc/w, the 
Strouhal number and v z (y^) the component of the velocity of the 
fluid normal to the wing at the seotion y^, the component by 
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assumption being constant along the entire section y^. Kussner, 
applying the acceleration potential of Frandtl also arrived at an 
Integral equation which resembles that of Cicala. ' If p equals 0, 
that Is, u(k) = 1, S(p|y - y^|) = 1 and both the equations of 
Cicala and Kussner became the equation of Prandtl which bolds for 
the case of the steady motion of the wing. lio solutions of the 
equations of Cicala and Bussner have as yet been found although to 
obtain these forms of the equations the authors Introduced Into the 
analysis many simplifying assumptions without which the equations 
would have been still more complicated. In Russia a paper on the 
theory of the wing of finite span In an unsteady flow has been 
published by M. Keldish but the latter’s paper contains the results 
of several stages of his work rather than a detailed presentation 
of his analysis. 

We see that notwithstanding the rather large number of papers 
devoted to the theory of the unsteady motion of an airfoil It must 
be confessed that In this domain the phenomena are still not 
completely Investigated. In addition to continuing purely 
theoretical Investigations, strictly scientific experiments are 
required both In order to explain the physical side of the phenomena 
and for the continued comparison of the results of theory and 
experiment. At the present time It is still far from evident that 
the physical picture which is assumed as the basis of the theory 
of the phenomena describes the latter completely. In particular, 
for example, it is entirely not clear how this physical picture 
must be changed for oscillations with finite amplitudes. Scientific 
experiments should aid in dealing with these problems. It Is 
desirable also to carry out a comparison of as many as possible 
different applications of stationary and nonstationary aerodynamics 
to the same phenomena. The theory of the unsteady motion of an 
airfoil Is a field In which many scientific conquests yet remain to 
be madg. 


Translation by S. Reiss, 
National Advisory Committee 
for Aeronautics. 
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Table | 


t 

at 

ac 

«c ow 

acf 

0 

0 

a 

a 

a 

1 

a 

1.105a 

1.649a 

2.718 a 

2 

2 a 

1 .221a 

2.718a 

7.389a 

3 

2a 

1.351a 

4.482a 

20.086a 

4 

4 a 

1.492a 

7.389 a 

54.598a 

5 

5 a 

1.649a 

12.183a 

148.413a 


Table 2 



s. 

M 2 

Hi 

c, 

n 

1 

•c"“* 

V, 

M/cen 

It 

Stork ......... 

1 " 

2.03 

0.123 

12.6 

11 

0.134 

P 1 9 eoa 

0.075 

0.67 

0.045 

31.4 

25 

0.056 

rook 

*0.130 

0.95 

0.0S8 

19.0 

12 

0.10B 

seagu II 

0.126 

0.96 

0.066 

13.9 

12 

0.076 
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Figs. 1,2 




Fig. I 




Fig. 2 
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